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Abstract 

The vibrating behaviour of thin structures is affected by prestress states. Hence, the 
effects of thermal prestress are important research subjects in view of ambient vibra- 
tion monitoring of civil structures. The interaction between prestress, geometrically 
non-linear behaviour, as well as damping and its coupling with the aforementioned 
phenomena has to be taken into account for a comprehensive understanding of the 
structural behaviour. Since the literature on this subject lacks a clear procedure 
to derive models of thin prestressed and damped structures from 3D continuum 
mechanics, this paper presents a new derivation of models for thin structures ac- 
counting for generic prestress, moderate rotations and viscous damping. Although 
inspired by classical approaches, the proposed procedure is quite different, because 
of (i) the definition of a modified Hu-Washizu (H-W) functional, accounting for 
stress constraints associated with Lagrange multipliers, in order to derive lower- 
dimensional models in a convenient way; (ii) an original definition of a (mechanical 
and thermal) strain measure and a rotation measure enabling one to identify the 
main terms in the strain energy and to derive a cascade of lower-dimensional models 
(hi) a new definition of "strain-rotation domains" providing a clear interpretation 
of the classical assumptions of "small perturbations" and "small strains and mod- 
erate rotations" ; (iv) the introduction of a pseudo-potential with stress constraints 
to account for viscous damping. The proposed procedure is applied to thin beams. 
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1 Introduction 



Ambie nt vibration monitoring (e.q. lWenzel and Pichlerl (120051 ). iBasseville et al. 
(120041 ) ) has now become widely accepted as an important tool for Structural 
Health Monitoring (SHM). But structural vibrations are affected by prestress 
states. In particul ar, thermal variations may cause very significant changes in a 
bridge spectrum (IPeeters and De Roecku200ll ; lFarrar et al.l . Il994j ) . Since tem- 
perature effects may be orders of magnitude larger than the effect of a damage, 
overlooking them prevents from a reliable damage detection based on vibration 
monitoring. Most attempts to eliminate thermal effec ts have favoured bl ind ap- 



proaches not taking advantage of predictive models (jSohn et al.l . 120031 ) . How- 



ever, successful endeavours to eliminate the temperature from subspace-based 
damage detection algorithms prove the relevance of relying on predicti ye ther- 
momechanical models yieldi ng the prestress state due to temperature (INasserl . 



20061 ; IBasseville et al.l . 120061 ). This suggests to deeper understand the way tem- 



perature interacts with structural dynamics and to revisit associated models. 
This paper steps forward in this direction. 

On the other hand, identification of the geometrically non-linea r behaviour of 
beam s or plates has recently g iv en rise to numerous pa pers (e.g. (IKerschen et al. 



20031 ). (lArgoul and Lei. 12003) . (IPerignon et al 1. l2003l)h often using the notion 
of non-linear modes (jRosenbergi . 119621 ) (jVakakisl . Il997l ). Non-linear dynamics 
of prestressed beams or plates und er going small strains an d moderate rota - 
tions is simulated in (lRibeirol . [200lh . (IPerignon et "all 120031 ). (lAmabilil . liooih . 
However, the dynamics used in these contributions seems to be based on some 
"historical" assumptions, whose justification is often skipped. 

The purpose of this paper is then to provide a new viewpoint on this classical 
subject, where a very large number of contributions, sometimes not clearly 
related, have been superposed over the years. We make clear the series of sim- 
plifying assumptions leading from the 3D continuous thermo-elasticity theory 
to the equations governing the dynamics of thin structures with prestress 
and thermal field, under the assumption of "small strains and moderate ro- 
tations". Attention is paid to the definition of the range of validity of these 
classical equations, introducing the notion of strain-rotation domains. More 
precisely: (i) a Hu-Washizu functional including suitable stress constraints 
and the associated Lagrange multipliers is defined: the constitutive law, the 
strain-displacement relationship as well as the dynamic equilibrium equations 
are then derived by imposing the stationarity conditions. Stress constraints 
characterizing thin body theories are introduced in the 3D models in view of a 



* Corresponding author. Tel: +33 1 64 15 37 80, Fax: +33 1 64 15 37 41 
^Acknowledgements: this work has been partially supported by the European RTN 
SMART SYSTEMS, while the first author was visiting the LCPC 
Email address: erlicher@lami.enpc.fr (Silvano Erlicher). 



2 



convenient derivation of ID beam models, (ii) The most important configu- 
rations characterizing the structural behaviour are clearly identified, (iii) The 
assumption that the prestressed configuration has a known geometry, ado pted 
as Lagrangian reference configuration (see e.g. ( Geradin and Rixen . 19951 )). is 
removed. When the pre-stress field and the original geometry of the structure 
are very simple, this hypothesis is convenient. However, for more general con- 
ditions, e.g. when the prestress state is time- dependent or related to a general 
thermal field, this assumption appears to be over-simplified. Hence, a differ- 
ent analysis, where prestressed and Lagrangian reference configurations are 
distinct , seems to be more suitable. The equilibrium equation governing the 
dynamics of a structure subjected to a generic prestress is defined as the dif- 
ference between the Lagrangian equations at the dynamic configuration and 
at the statically prestressed configuration, (iv) A rigorous formalization of 
the assumption of "small strain and moderate rotations" is provided. To this 
end, measures of the strain amplitude (symmetric part of the displacement 
gradient) and of the rotations (skew-symmetric part) are introduced in the 
3D continuous framework, where as standard a pproaches use the aspect ratio 
as the governing parameter (e.g. ( Ciarlet . 1980h ) . The new solution-dependent 
measures enables one to define strain-rotation domains in which the leading 
terms of the strain energy are clearly identified. 



Damping, due to internal friction or other dissipative phenomena, should care- 
fully be taken into account in the structural analysis. Several damping models 
exist, like general linear damping or viscous proportio nal or non-proportional 
linear damping; see e.g. ( lAdhikari Woodhousd . 120011 ). However, the link be- 
tween these structural damping models and the corresponding dissipative ma- 
terial behaviour, described by a given strain-stress law, does not seem to be 
clearly established in the literature. Here the beam dynamic equations with 
linear viscous damping are derived for the case of "small strains and moderate 
rotations", thus filling the gap. 



The different configurations characterizing a vibrating structure subjected to 
static and dynamic loads and to a thermal field are defined in Section 2. 
Then, the thermo-elastic constitutive rule is presented in Section 3. Physi- 
cal linearization is considered. Section 4 introduces a modified Hu-Washizu 
functional accounting for stress constraints. The corresponding stationarity 
conditions are discussed. In Section 5, two global measures for the strains and 
the rotations are introduced and used to suggest approximated expressions of 
the strain energy, each approximation being valid in a strain-rotation domain. 
2D finite element simulations substantiate the approximation of the strain 
energy. Section 6 introduces a dissipative stress tensor in view of modelling 
damping effect. Section 7 leads to a general expression of the weak equilib- 
rium of a continuum subjected to a general static prestress. In Section 8 the 
previous general procedure is used to derive the Euler-Bernoulli beam equa- 
tions for small strains and moderate rotations in the undamped case, while 
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the damped case is treated in Section 9. After the Conclusions, the Appendix 
explains how to compute the Lagrange multipliers. 



2 Configurations of a structure 



The following configurations (see Figure [T] and Tabled]) may be distinguished: 

(1) Voo is the geometric reference configuration, where the displacements are 
assumed to vanish. For a beam, Voo is the straight configuration. Let X 
be the reference position of a material point. 

(2) Voi is the initial configuration of the body, where the temperature field 
Toi is constant and where no external force is applied. This configuration 
is relevant in the case of geometric imperfections, where V^ differs from 
Voo- In this paper, however, we assume Voi = Voo, Uoi = x j — X = 0. 
The stress may not vanish since a self-equilibrated stress n j may exist. 

(3) Vq is the equilibrium configuration under static external loads f , go and 
a temperature field T resulting in a displacement U , a relative displace- 
ment between Voi and Vq equal to Uo — Uoi and a stress field n — n j. 

(4) V± is the instantaneous configuration of a vibrating structure subjected to 
the static loads f , go, to the temperature field T and to dynamic volume 
and surface forces fi and gi. The relative displacement field between Vo 
and V\ is indicated by U i = Ui — U = x x — x . In practical situations 
fo , go and To may vary with time, e.g. on a daily period, but their 
variations are supposed to be very slow with respect to free vibrations of 
the structure. This decomposition of a given external force into a static 
part and dynamic part is convenient but not unique, e.g. it may depend on 
the time interval considered if two different static loads are successively 
imposed. Therefore, the total response of the structure does not split 
uniquely in a static component and a dynamic response. 



3 Constitutive law and physical linearization 



Based the objectivity principle, (see e. q. lMandell (11 9 661 . p. 602)), the Helmholtz 
free energy \1/ can be defined as a function of the Green-Lagrange strain tensor 



E(U) 



<9X " 1 " \dX J 



T 2 \dX) 



■(1) 



V s (U) + ±V (U) T ■ V(U)(1) 



and of the temperature T, viz. \1/ = \& (E, T). As usual, • indicates the dot prod- 
uct and U is the displacement measured with respect to the given reference 
configuration Voo- When ||E|| <C 1 and the temperature variations are small, 
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^ (E,T) can be approximated by a truncated series expansion. Truncating at 
the second order around the given initial configuration Voo, characterized by 
E = and T = T 0i (physical linearization) leads to 

(E,T) = n 0i : E+±E : D : E - E : A (T - T 0i ) 

2 

—Soi (T — T 0i ) — (T — T 0i ) 



where : denotes the doubly contracted inner product, D the fourth order tensor 
of the elastic constants, T — T 0i the temperature variation and A a diagonal 
second order tensor accounting for thermal expansion, Soi the initial volume 
entropy and c £ the specific volume heat [J m~^K~ l \. When a St. Venant- 
Kirchhoff material is considered, one has 

D =A1 ® 1+2//I and A =q (3A + 2/i) 1 (3) 



where (g) is the outer tensor product; 1 is the second order identity tensor; I 
is the fourth order identity tensor; A and \i are the Lame coefficients, a is the 
thermal dilation coefficient. As it is well-known, the following identities hold: 

= — = - (d) 

(! + !/)(! -2i/)' M ~ 2(1 + 1/)' 1 j 



where E is the Young modulus, v is the Poisson ratio. In this case, the 
Helmholtz energy (j2J) becomes 

* (E,T) = n 0i : E+|A (tr (E)) 2 +juE: E— a (3A + 2fi) (T - T 0l ) tr (E) 

(5) 



— Soi (T — T 0i ) — (T — Toi) 2 



Eq. ([2]) leads to a constitutive law linear with respect to E and T — T 0i : 
n = || = n" d = n 0l + D : E - A (T - T 0l ) (6) 



where II is the (second) Piola-Kirchhoff symmetric stress tensor, the index nd 
means non-dissipative. Eq. ffih also reads 



n = U nd = U 0i + Xtr (E) 1+2// (E) - a (3A + 2/i) (T - T 0i ) 1 (7) 



For isotropic materials having elastic non-linear constitutive behaviour, a di- 
rect generalization of (JTj) can be defined. It is called second order elasticity 
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( jMandell . Il966l . p. 607), supplemented here by the self-equilibrated stress Ylof. 



n = n 



nd 



n 0i + 



A tr (E) + |tr (E 2 ) + W (tr (E)) 2 - a (3A + 2/i) (T - T { 



On 



+ 



A'tr (E) (T - T 0i ) + a' (T - T 0i ) 2 1 + 2/i E + A tr (E) E + 2/i'E (T - T w ) + C E • E 



where A, B, C, A', //' and a' are the material parameters introduced in addition 
to the usual ones A, /i and a. 



Remark: the temperature variations can be considered small when 



D : A (T — T c 



Oh 



< 1 



(9) 



viz. when strains associated with them are small with respect to unity, simi- 
larly to strains fulfilling condition ||E|| <C 1. For isotropic materials, the tem- 
perature variations are small if \a (T — T 0i )\ <^ 1. This condition is consistent 
with the assumption of physical linearization. 



4 A Hu-Washizu functional with additional stress constraints 



In this section, a special representation of the problem at hand is introduced, in 
view of taking into account stress constraints in a three-dimensional framework 
without forgetting about compatibility issues. A Hu-Washizu (H-W) func- 
tional depending on the stress II, the strain measure E and the displacement 
U, considered independent is supplemented with constraints on the stress, in 
order to model thin bodies. Let U , T and F denote spaces of smooth enough 
displacement fields, of symmetric second order tensor fields and scalar fields 
of Lagrange multipliers, respectively. The proposed H-W functional reads 

J ff _ w (ir,E*,U*,A*) =/ Voo *(E*,T) dV-f Voo II:(E*-E(U*))dV 

- J Voo f • UW - ! dVo0 a g • V*dA - J 9Vm u ([(1 + V (U*)) -IP] -N) • ( U* - U) dA (10) 

- Ivoo nti KR k ■■ n w 



with IT G T, E* G T, U* € U and A* G F" A , where n A is the number of 
stress constraints and tj is the given displacement on the boundary dV 0QjU = 
dV 00 — <9Voo, CT . Note that Jh-w is a functional defined on 3 + ri\ fields. At a 
solution (^11, E, U, Aj of the static problem, the functional above satisfies the 
stationarity with respect to the 3 + n\ fields. The last term leads to linear 
stress constraints: 

R fc :n:=(N fc ®N fc ):n = (n-N fc )-N fc = (11) 
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where is a vector and R& is a second order constant tensor. Eq. (TTTT) im- 
posed on II carries over to an equivalent condition on the Cauchy stress cr. 
The use of this kind of constraints to derive beam equations from 3D elas- 
ticity is illustrated in Section 7. The case of a general 3D structure with 
no constraint can be addressed by formally setting = and R& = 0. 
The well-known impossibility to derive beam or plate equations from purely 
kinematical assumptions in the three-dimensional equations of elasticity in 
pure displacement motivate the introduc ti on of internal constraints as in 
( INardinocchi and P. Podio Guiduglil . Il994j ). (ILembo and Podio Guiduglil . 120011 ) . 
These constraints enable one to use a pure displacement approach while mim- 
icking the averaging process underlying the convergence of the equations of 
elasticity when the aspect ratio tends to zero. The derivation of lower-dimension al 
mode ls without assumptio ns can rely o n T-convergence (I Acerb i Buttazzol . 



1986), asymptotic analysis (jCiarletl . Il980l ). or energy methods (IBabuska et al. 



19921 ). This paper aims at deriving the equations governing the evolution of 



thin structures subject to prestress states without above mathematical appa- 
ratus. 



4-1 Action functional and stationarity conditions 



While the stationarity of the H-W functional leads to statics, the action 
functional leads to elastodynamics over a time interval [0, tf\. Let us set 
V = {U (X, t) : U (-, t) G U Vt G [0, t f ]}. The kinetic energy T (u*) = \ J Vw p 00 U* ■ UW, 
combined with the H-W functional enables one to define the action 

v (it, e*, u*, a*) = J (r (u*) -j h-w (n*, e*, u*, a*)) dt (12) 





Hamilton's principle is given by 8T> = 0, where small variations (SH, 5E,5U ,5 A) G 
T x T x V x F nA are considered. The stationarity operat or 8 (•) is intended as 



isochr onous, viz. 8t = and 8 J f Cdt = J f SCdt (see also lQuadrelli and Atluri 



(119991 )). Hamilton's principle leads to the stationarity conditions 



I. 
2. 

3. 



V7c = 1, Rfc : II = in Voo 

E = E (U) - Y2t\ AfcRfe in^oo 
U= U on dV n u 



(13) 



n 



9E* 



in Vr 



00 
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In the physically linear case, Eq. ([131-3) corresponds to dHj). If isotropy is as- 
sumed, Eq. (T13I-3) is equal to (|7j). The first equation expresses the stress con- 
straints, the second one defines the strain E, which differs from E (U) due 
to the stress constraints. Moreover, II depends on the strain E accounting for 
the Lagrange multipliers and not on the standard strain measure E (U). The 
stationarity condition of the action T> with respect to displacements leads to 

V<5U G V, /*' \j Voo n : SE (U,5U) dV - J Voo f-SUdV - J dVo0 a g-SVdA 



Io f \ldv 00 , u ([(1 + V (U)) -H] -N) -6UdA + J Voo pooU-aW 



dt = 



dt 

(14) 



where it has been used the identity U= U on <9Voo,u, given in (TIB"]) . Moreover, 
the virtual strain is defined by 

5E CU,5V) = ■ 5U 

v ' ' du u*=u 

= V s (5U) + | V (5U) T • V (U) + |V (U) T ■ V (SU) 



Integrating by parts in time leads to the virtual works principle at every t 
W t (n, U, 5V) + W e (f , g,5U) = W a (V,5U) V5U G V (16) 

where 

(n, U, 5U) := - J Voo n : 5E (U,5U) dV, W a (U,5U) := J Voo Poo tj-5UdV 
W e (f , g,5U) := f Voo {-5UdV + f 9Vo0a g-5UdA + J dVwu ([(1 + V (U)) n] -N) -5XJdA 



denote the virtual work of internal, inertia and external forces, respectively 
Finally, by integration by parts, one obtains the corresponding strong form 
equation and the boundary conditions on dVoo,a- 



div 1+V U -n +f =PooU inV^oo 
4. <^ (18) 
;i + V(U))-n-N = g on^oo,a 



with the initial conditions U (t = 0) = U and U (t = 0) = V in Vqo- 



(17) 
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5 Physical linearization, strain-rotation domains and dominant terms 
in the strain energy 

Let U, E, II and denote the fields defining the equilibrium at time t. The 
strain energy density associated with the strain field E = E (t) reads 

F(e,t) = J #(E,T) dV (19) 

Voo 

where \1/ (E, is given by (J5]). The displacement gradient writes 

V (U) = V s ( U) + V sk (U) := e (U) +r (U) = e + r (20) 



where the symmetric tensor e is called small strain tensor, while the skew- 
symmetric part r is called the small rotation tensor. For brevity, hereinafter 
they will be referred to as strain tensor and rotation tensor. In the same way, 
the strain increment ffl5|) writes 

5E(U,5U)= V s (<JU) + (<5U) T • (e + r) + ~ (e + r) T • V (5U) (21) 



5.1 A global measure for strains and for rotations and a hierarchy for the 
strain energy terms 



In view of finding out the dominant terms in (fT9]) . let us introduce the solution- 
dependent scalar parameters r\ and p as follows: 

V = Ve + VAT, Ve = 7= \\e || i2 (y oo) , ^AT = ^ ||« (T - T 0l ) l^y^ 



^ ''~ V^VWi H r lli2(Vbo) ~~ \Voo Ivoo 2 ri fv dV 



where ||'|Il 2 (Voo) * s ^ e -^2" n orm associated with the domain Voo- The strain 
measure r\ is supposed to be small. Observe that: (i) p is only defined when- 
ever 77 > 0. Thus infinitesimal rigid body motions at initial temperature are 
excluded from the following analysis; (ii) If the rotations are not large, r is a 
rotation tensor and the quantity rj p is a global measure of the rotations. Hence, 
p can be interpreted as a global index of the relative amplitude of rotations 
with respect to strains. As a result (see Eq. (Tj[|)) 

E (U) = O (17) + O (V p ) + O (v 1+P ) + O (V) (23) 
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Under dynamic conditions and provided that the vibration frequency is bounded 



one also has 



dt 



O (77) and 



£ 2 (Vbo) 



O (r) p ). In addition, let us 



make a consistent hypothesis on the self-equilibrated prestress n j: 



n r 



nj? + ng p) + ng) + n£ +p) = 0(^ + fo*) + o (^) + o tf) (24) 



Since the Lagrange multipliers A& fulfill the relationship for all k — 1,71a 



Rt- : II — Rj. : 



"A 



n , + D: E(U)-^A ; R -A (T - T«) 



(25) 



then, provided the stress constraints R& are linearly independent, A= (A fc ) , 
fc = l,riA> is the solution of a linear Gram system and splits as follows (see 
Appendix): 

A=A (1) +A (2p) +A (1+p) +A (2) =0(t]) + (rf p ) + O (V +p ) + O (rf) (26) 



where A*- 1 -*, A^ 2p \ A^ 1+p ^ and A 1 - 2 -* can be evaluated analytically in simple 
cases. This implies that the Lagrange multipliers, E (Eq. (TT51-2)) and E (U) 
split in terms of the same order. 



From Eqs. (1131) and (I2"6"j) . the stress splits in four terms: 



(27) 



where 



n« = n« + D : (e - Er= A i aJ 1} R) -A (T - T w ) = O ( V ) 

n (2p) = n (2p) + J) . l r T r _ j) . A (2P) R;=0 (^2p) 

n (i+p) = n g+P) + D : (|r r ■ e + \e T ■ r) - D : £?4 A{ 1+p) R = O (r7 1+p ) 
n( 2 ) = ng> + D : \e T ■ e - D : £^ Ap } R=0 (r/ 2 ) 



(2* 



Finally, assume that 

0<r/<l p > 



(29) 



Since 77 < 1 , the terms with the smallest exponent are the largest ones. The 
comparison between different terms must be done for a given p- value. If p > 1, 
viz. when r is of the same order of magnitude as, or smaller than, e (case of 
small rotations), all non-linear terms in the Green-Lagrange strain E and 
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in the Lagrange multipliers can be discarded and the standard condition 
of small transformation is retrieved. The exponent p < 1 is introduced in 
order to account for the possibility of having rotations larger than strains. 
The complete expression of the strain energy reads (recall Eq. (121), where the 
terms that depends only on the temperature are omitted for brevity): 

T (E, t) = f Voo v& (E, t) dV = f Voo (n 0l — A (T — T 0l ) + |D : e) : EdV 
= J VoQ (n 0l — A (T — T oi ) + |D : (E- E^i A,R,)) : ( E- A,R,) dF 



i.e. ^ (E, T) = ^ (E, T) (E, T) , where 

.F (E, T) = J yoo (n 0l — A (T — T 0i ) + |D : E) : Ed^ 

JF A ( E , T) = — J Voo (n 0i — A (T — T ,) + |D : E) : (£J* A Z R,) dV 

~ Iv 00 |D : (£74 W : Erfy + |D : (E^i A,R) : (£& A/R;) 



More precisely, 



F (E, T) = J Voo ( ng> — A (T — T 0l ) + |D :e ) : (e +\r T ■ r) dV 



+ ivbo 
+ JVbo 
+ JV00 
+ JVbo 
+ ivbo 



n 



(i) 

Oi 

(2P) , 
Oi "T 

(1+P) 
Oi 

r(l+P) 
l 0i 



A(T-T 0j ) + iD :e):(|r r 



e + \e T - r- 



-l£ T 

2 C 



ID: 



ID: ir 



IT . 
2 

1„T 



e +|r T • r- 



2 r 



£ + l£ T ■ r- 



dV 



1 T 

■5 e ' £ 



+ id 



2 

1 T 

-r 
2 1 



£ +l£^ 
£ + l£ T ■ 1 



Ilg> + ID: (l£^. £)):(£ +Ir^.r- 



+Ir T ■ r 

2 e ~ 2 



Ir^.£ +i£^ 



d\/ 



dV 



(30) 



r+l£ T ■ £ 



d^ 

dy 



Eq. (1301) shows all terms of the strain energy not depending on A. The expo- 
nents Si of these energy terms are indicated in Figure [2] and in Table [2] The 
strain and rotation measures 77 and rj p lead to a hierarchy between the different 
terms. For every p— value, the smallest exponents are si = 2 , s 2 = 1 + 2p, 
S3 = 4p (see Fig. [2]). Then, the dominant terms are 

(E, T) = J Voq (ng> — A (T — T ,) + ID :£ ) :£ dV = O (r? 2 ) 
^(i + 2 P ) (E) T) = /voo ( n (i) _ A (T — T 0J ) + ID :£ ) : ■ r) dF 

+ /voo K P } + |D : (Ir* • r)) :s dF = O ( V 1+2p ) 
( E , T) = / yoo (n^ + ID : (Ip^ • r)) : ■ r) dV = O {rf*) 
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Likewise, the dominant contribution to Tk splits in three terms: 

(E, T) = - J Voo (ng) — A (T — To,) + |D :e ) : E^i A^R^V 

- Jy 00 I ( D : Sr= A i a[ 1} R) :e dV + J Voo \ (D : Efii A^R,) : Ya=\ A, W R|dV = O 

4 1+2p) (E, T) = — (n« — A (T — To,) + |D :e ) : Af^dV 

- K P) + |D :|r^ ■ r) : E?ii A«R^ 

- /voo I ( D : A, (1) R) ^ • rdV - J Vqo § (d : A^B,) : edV 
+ §2 (D : EFii A^R,) : E™=i Af p) R^y = O ( V 1+2p ) 

F<*» (E, T) = - / yoo (l4*> + |D :|r^ • r) : Efii A^R,dV 

-J Voo |(D:EriiA; 2p) RO:|^.rdy 

+ I Voo \ (D : E"=i Af P) R) : E™=i Af^dV = O (r/ 4p ) 

In summary, one can write 

T (E, T) = ^ (E, T) + (E, T) + J2 (V Si ) (32) 



with 

jib) ( E> T ) = jr(2) ( E) T ) + jr(i+2 P ) ( E) T ) + jr(4 P ) ( E? T ) 
(E, T) = j| 2) (E, T) + j| 1+2p) (E, T) + j| 4p) (E, T) 



A strain energy containing only the three dominant terms, viz. JF = JF^ (E, T) + 
^ r | fe ' ) (E,T), can be directly derived starting from the Hu-Washizu functional 
(HO I) , provided that the strain measure is chosen as follows: 

E (U) - EW (U) = e (U) + X (U) , X (U) = ±r (U) T ■ r (U) (33) 
The virtual works of the internal forces (see (JIT]) . (j2*Tj) and (1281) ) read 

-nf } (n, u, 5U) = j Voo n«:V s (*u) dv 

-w\ 1+2p) (n, u, su) = j Voo n^:Sx (u,*u) dv + j Voo n^ : v s (su) dv (34) 
-wf 4p) (n, u, su) = j Voo n^:5 X (u,*u) dv 
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where 5\ (U,5U) = \ (r (5U) T • r (U) + r (U) T • r (5U)) . 



5.# Definition of "strain-rotation domains" and of approximated expressions 
of the strain energy 

Eq. (1321) reveals the dominant energy terms. However, it is also interesting to 
find the strain-rotation domains (depending on 77 and p) where the approx- 
imation of the exact strain energy by a simplified expression is acceptable. 
Three domains are of special interest: 

(a) The domain M a where 

T (E, t) ~ (E, T) + (E, T) = (rj 1 ) (35) 
Eq. (135]) represents the situation usually called of small perturbations. The 

(2) 

corresponding virtual work of internal forces reduces to Wj (II, U, SU); see 
Eq. ( 1341) . As it is seen hereafter, a suitable name for H a is domain of small 
strains and relatively small squared rotations (the ratio r/ 2p_1 is small). 

(b) The domain H& where 

T (E, T) ~ (E, T) + J-f (E, T) = O (if) + O (r/ 1+2p ) + O (r] Ap ) (36) 

ie. three energy terms are considered together, of order 2, 1 + 2p and 4p, 
after the discussion of the previous Section. The virtual work reduces to 

wf } (n, u, sv) = wi 2) (n, u, <ju) + wi 1+2p) (n, u, sv) + wi 4p) (n, u, su) 

This region is the domain of small strains and relatively moderate squared 
rotations (the ratio rf v ~ x — 1), called moderate rotation domain for brevity. 

(c) The domain H c where 

T (E, T) ~ ^ (4p) (E, T) + Tf v) (E, T) = (r] 4p ) (37) 

and Wi (n, U, SU) -> (n, U, 5U) = O (r] 4p ) , ie. the energy term of 

order 4p is larger than all the others (domain of small strains and relatively 
large squared rotations, i.e. the ratio rf p ~ l is large). 

In order to draw these domains, let us introduce a small positive number 
( C 1, for instance ( = 0.01. Then, the situations p < | and p > \ can be 
distinguished, since for p < ~ the dominant term is of order 4p and in the 
other case the term of order 2 is the largest (Figure [2]). When p > |, look for 
the conditions on 77 and p such that all energy terms are small compared with 
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the second order term. The following system of inequalities defines M a : 

U a : ^-<C foralH = 2,n (38) 

Every condition enjoys a simple log- log representation, (x = logr] p ,y = log rj), 
since r] Si ~ 2 = r] ClP+di — > log77 Si ~ 2 = log7? CiP+dl = CjX + diy < log£, where c, 
and dj are integer numbers. For instance, when Sj = s<i = 1 + 2p < log(, one 
has 2x — y < log£, which relates ?/ = log 77 with 2 = logr] p . This condition is 
associated with a line which is the bottom limit of the domain H a ; see Figure 
If Si = 3, one has q = 0, di — 1, leading to the condition y = log// < log£, 
associated with an horizontal line as shown in same Figure. This means that 
the ratio between the energy term of order 3 and that of order 2 is smaller 
or equal to (, provided that r\ < £. Still assuming p > | , one can define 
the conditions such that all the energy terms are small compared with that of 
order 2, except those of order 1 + 2p and 4p: 

H' 6 :^<C Vz = 4,n, C< J V:=Pi<l, Cg^y :=P2<1 (39) 



The third relationship shows that inside this domain, the ratio between the 
term of order 4p and that of order 2 may be either small, equal or greater 
than (. When p < ~, the conditions such that all the energy terms are small 
compared with that of order 4p read: 

H c : ^ < C for all i = 1,2 and 4, n (40) 



V 



This set is associated with the approximated energy (1371) . If the terms of order 
1 + 2p and 2 are not small, one has 

K ■■ \ < C Vi = 4,n , C < %z- = - < 1, ti\ = -<l (41) 



The set H& = W b U H^' is associated to the approximated energy fl36l) . Inside 
this domain, three energy terms (of order 2, 1 + 2p and 4p) are retained. The 
ratios p\ and p 2 and their inverses (Eqs. (139]) . (jHJ) give an estimation of the 
relative amplitude of these dominant terms. 

Observe that H c is extended to large rotations (77 4p close to 1) because we have 
compared the energy terms of the first four rows of Table [2] characterizing the 
physically linear and isotropic constitutive law (J7J). However, it is of interest 
here to determine the conditions under which the physical linearization is 
a reasonable approximation of any real material behaviour. It is expected 
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that this is the case when the rotations are not too large. In order to derive 
sharp bounds on the rotations, the quadratic constitutive law (JHJ) has to be 
considered, and the associated strain energy must be computed. This leads to 
(see also (127]) ) 



n (phys.non-lin.) = jjiphys.lin.) + ^ + q ^2p) + Q + Q 

T = (t] 2 ) + O (t] 1+2p ) + O (t] 4p ) + E?=t O (rf <) 



i.e. some energy terms have to be added to those of the physically linear case, as 
indicated in Table [2j Then, the same conditions as in (!38l) - (|4T!) are considered, 
with n substituted by n n i in , in order to compute the conditions under which all 
the physically non-linear terms are small compared with the dominant ones, 
either of order 2 or 4p. These conditions define the regions depicted in Figure 
[3b: they are the strain-rotation domains where the physical linearization is 
admissible. An important difference with respect to Figure [3^ is that H& and 
H c are bounded by the vertical line x < | logC, equivalent to rj 2p < (, viz. the 
squared rotations, not only the strains, must be small. It can be proven that 
this limitation derives from the condition ^ < (, imposing that the term of 
order Qp associated with the physically non-linear law remains small compared 
with the term of order 4p, dominant for p < 1/2. 

The condition of having small Green-Lagrange strain reads: ||E|| = O (rj) + 
O (r] 2p )+0 (r] 1+p ) + (rf) < £ <C 1. It easy to identify in Figure[3]the domains 
in the strain- rotation plane where the first, the second and the fourth term of 
E are less or equal to (. It can be also proven that the condition of having a 
small third term, i.e. i] 1+p < (, is fulfilled in the three sets H tt , and H c . 

Other approximations retaining at least four energy terms are possible. How- 
ever, (a), (b) and (c) define situations often discussed in the literature and 
for this reason the present analysis is restricted to them. Case (b), collecting 
three terms, is formally more complex than the others and is discussed in 
detail hereinafter. 



5.3 Numerical example: "exact" values of rj and p and comparison of the 
exact and "simplified" energy maps 



Consider a problem of plain stress elasticity (n^ m , the dimension of the prob- 
lem, is equal to 2). The corresponding conditions on the Piola-Kirchhoff tensor 
are Tlxz = ^-yz = ^zz = 0. A St.Venant- Kirchhoff material is chosen with 
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Iloi = and T = Toj. Then 

* = * (E) = |E : (Al ® 1+2/il) : E 
II = II (E) = 2/i E+A tr (E) 1 = 



IJxx Iljfy 



(43) 



where A = 7^ A due to plain stress assumption and 



E(U) 



du 1 ^ 9u | 9v 



ax 

3u I dv \ 
dY ^~ dX ) 



2\dY 



dX 



dv 
dY 



du du I dv dv du du i i9v 3v 

ax ax ax ax ax ay ^ ax ay 

au 8u I dv dv du du i 3v 3v 

ax ay ax ay ay ay ay ay 



is the 2D Green-Lagrange strain with U = [u, v] T . The relevant H-W func- 
tional is given by (flOj) . without the term depending on A* k and with \1/ given 
by (H31-1). For a static problem with the external volume force f = [fx,o, /y,o] T 



and surface force g 
reads: 



[gx,o, 9y,o] j the weak form of the equilibrium equation 



R(U,5U) = -/ Vbo n(E(U)):5E(U,5U) dV 

+ J Voo f -5U dA + J dVo0a g -5 U dA = for all 6U GV 



(45) 



with SU = [Su, Sv] T and SE (U,8U) defined by the right-hand side of Eq. (fT5|) . 
Eq. (145p is a non-linear partial differential equation, which can be discretized 
by a standard finite element method. The s tandard Newton alg orithm has 
been implemented in the code FreeFEM++ (IDanaila et all 120031 ): given the 
displacement U n at iteration n, the increment w is computed by 

w eU such that R (U n , 5U) + SR (U n , SU, w) ~ for all 5U GV 



with 



5R(U„,5U,w) 



- / n n (SE (U n , w)) : 5E (U n ,5U) rffi 
/ n n (E (U n )) : 6 2 E (SU, w) dfi 



5 2 E (5U, w) = |V (5U) J ■ V (w) + ±V (w) J ■ V (5U) 



Set U n+ i = U n + w. Repeat until MWml j s sma ll enough. 

H u "lli 2 (Voo) 

The structure examined in this example is a parallelepiped beam lying in the 
XV-plane, the dimensions are b = h = lcm and Lqo = 50cm. The beam 
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is clamped at both ends and the volume load is [fx,o, /y,q] = [0, — |/y,o|] 
daN/cm 3 . Exploiting the symmetry of the problem, only a half-beam is meshed, 
with the boundary conditions [u = 0, v = 0] for Y G [—h/2, h/2] and X = 
and u = for Y G [— h/2, h/2] and X = L 00 /2. The material parameters 
(steel) read E = 2100000 daN/cm 2 , u = 0.28, // = 820312.5 daN/cm 2 , 
A = 638020.8 daN/cm 2 . A mesh of triangular elements has been chosen, with 
two elements inside every cell of a regular grid of 15x375 squares. The finite el- 
ement space is of PI type. The numerical simulations give the results collected 
in Tables [3] and H, where v max = v (X = L 00 /2), x = log 10 (rj p ) , y = log 10 (77) , 

V = y , Pi = V' P 1 ^ aIld 



^ P y/Ww) ll r lli 2 (Vbo) V bhL 00 /2 




(46) 




The last column of Table H] provides a global estimation of the energy er- 
ror between the exact and the approximated energies T = J (E) dV and 
dV. Observe the maximum absolute value of the strain e TT , 
reported in the third column of Table [3j in the first four cases it is less than 
0.002, which is the limit elastic strain for a steel having yielding stress ap- 
proximately equal to 4200 daN/cm 2 = 420 MP a. For these situations, the 
material is truly physically linear. Conversely, when /y = —6 daN/cm? , see 
the last row of Tabled the maximum absolute value of e xx is larger than 0.002. 
Hence, the physical linearity is truly fulfilled only for steels having a greater 
yielding stress. The x — y coordinates of the strain-rotation points associated 
with each fyo value are reported in the seventh and eighth columns of Table 
[31 The corresponding graphical representation is given in Figure [3]-b, depicted 
assuming ( = 0.01. The map \& (E) , viz. the strain energy density at the 
static equilibrium Uo is shown in Figure HI The energy density corresponding 
to the approximated strain energy \1/ ^E^) is given in Figure [51 The relative 
difference of energy density between the two cases is illustrated in Figure [6j 
With a surface load g = [gx,o,9Y,o\ = [0, — |fly,o|] daN/cm 2 and f = 0, the 
results are similar, as one can see from Tables E] and Figure [TJ This confirms 
that 77 and p are not too sensitive to the load distribution. 

6 Dissipative stress, Hu-Washizu functional and damping pseudo- 
potential with stress constraints 

The linear law can be generalized by adding to U nd a dissipative term: 
n = U nd + U d (47) 
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The index d indicates the dissipative part of the stress. Let 

- Ivoo nrf * : (f * -§)dv- f Voo Efe=i KR k : n d *dv 



be the functional associated with a dissipative stress. It is assumed that it 



depends on the dissipative stress U d G T, the generic strain flow G T 
and A£ G F, i.e. the Lagrange multipliers associated with the constraints im- 
posed on U d . Moreover, the actual strain flow plays the role of additional 
parameter: for this reason it is separated from the main variables by the semi- 
colon " ;" , instead of the comma. The actual strain flow is computed from the 
problem associated with the Hu-Washizu functional defined in Eq. (1541) . The 
scalar non-negative and convex function is called pseudo-potential or dissi- 
pation potential. A classical definition is 4> = : F : ^ , i.e. a quadratic 
function. For an isotropic material, one has F =A^1 <S> 1+2/^1 , where A^ and 
fid are analogous to the Lame constants A and fi. The stationarity conditions 
imposed on fT4"51) lead to following strong form expressions 



a. Wk = 1, n A R fc : IT = in Vq 



oo 



b - IF ~ W ~ Sfc=i AfcRfc in Voo ^ 



n 



9^ 



F : § in V 00 



dE _dE 

dt dt 



The first equation indicates the stress constraints imposed on n^, the second 
one shows that the strain flow governing the dissipative behaviour is not equal 
to the time derivative of the strain when A^ ^ 0. Finally, the third equation 
is the constitutive law for the dissipative stress, obtained from the pseudo- 
potential 0. The constraints on the dissipative stress read 



Rfc : II — Rz 







(50) 



where Aj, I = 1, n\ are known from the analysis of the non-dissipative part of 
the stress. Following the procedure indicated in the Appendix, one can prove 
that A = (Afc) , k = l,n\ is the solution of a linear Gram system. Therefore 

\=\v + \(2p) + \a+p) + \w = 0(^ + ( v 2p ) + o ( v 1+p ) + o (V) (5i) 

Using Eqs. (14"§I) and (15T|) . one obtains 

YI d = YI d ^ + n "- 2 ^ + lT^ 1+p ) + Il d ^ (52) 
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where the following four terms of different orders are distinguished: 



n d(1) = F : (f - Efe=i ^R, - Efc=i a£°b*) 

n ^) = F:(|^.s +|r^.^ + |^.f + |^.r) (53) 

- F : (e& + E3b=i Al 1+p) R fe ) 



In order to define the non-dissipative part of the stress, as well as the equilib- 
rium equation, the following Hu-Washizu type functional is introduced: 



J H - W (lF d *, E*, U*,Al; n d ) = J Voo * ( E *, T) dV - J Voo II nd *: (e* - E (U*)) dV 

+ n< " e (u*) rfv - / yoo f • uw - g • u*^ 



IdVoo, 
IV 00 



;i + v (u*)) • (n nd * + n d )] -n) • (u* - u) 



(54) 



Eq. (1541) should be compared with (TlOT) . An attentive reader can see that 
the functional depends on the non-dissipative part of the stress II nd G T 
instead of on the total stress II*. Moreover, an additional dependence on U d 
is introduced, where II d is the stationary solution of (jUJ). Stationarity imposed 
on Jh-w leads to the strong form expressions: 



I. 
2. 

3. 



R fc : U nd = , fc = l,n A 
E = E (U) - Efe=i A fe Rfc 
U=U 



n 



d*(E*.T) 



<9E* 



in V^o 
in V 00 

on dVoo,u 

in Vqo 



(55) 



E*=E 



Moreover, recalling ( 1T21) and imposing stationarity in the displacements, one 
obtains the following weak form equilibrium equation: 



(ll nd , U,5U) + W, (n d , U,5U) + W e (f , g,5U) = W a (U,5U) , V<5U G V 
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where W t (n nd , U,5U*) is defined in (ED, W e (f , g,5U) is also given in flU}, 
with n = U nd + n d and 

Wf = W< (n d , U,<5U) = - J U d :5E (U,5U) (56) 

Vbo 

The virtual work of inertia forces is the same as in the previous non-dissipative 
case. The dynamics of the system is ruled by (118I) . where II is given by (147j) . 
The case of moderate rotations is obtained by just substituting E(U*) with 
(U*) in Eqs. (|48p and fl54l) . This corresponds to the substitutions W; — > 
Wf } and Wf -> Wf b) , where (see Eq. (JMD) 



7 Accounting for a static prestress 

As already discussed, when a static prestress due external mechanical and/or 
thermal loading occurs, the structure passes from the state Vqq = Vqi to a 
state Vq. It is interesting to write the equations governing the equilibrium 
at the generic configuration V\ as a function of the unknown displacement 
Uoi = Ui — U , expressing the motion with respect to Vq, as illustrated 
in Figure [TJ This can be easily done subtracting the equilibrium equations 
established in the previous sections, written at V\ and at Vq. Both equilibrium 
conditions at Vq and V% should be written, together with the other expressions 
coming from the stationarity of the relevant H-W type functional. For the sake 
of simplicity, only the weak form of the dynamic equilibrium is reported in the 
analysis of this section. As seen above, the virtual work of the internal forces 
is indicated by W, in the general case, and by wj^ in the case of moderate 
rotations. All the equations of this section are written as function of Wj, and 
then refer to the general case. However, the formal substitution of VVj at the 
place of Wj gives the equations for the moderate rotation case. 

The static problem defining the prestressed configuration Vq reads 
Find U GU such that for all 5XJ GV 

< w, (n , u ,5U) + w e (f , go , su) = o , n =n nd (u , a m ) (58) 

U = tJ on dVoo 7U 

where A^^ are the Lagrange multipliers associated with the static problem and 
the prestress n accounts for the temperature field T . The dynamic problem 
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denning the generic configuration V\ reads 



Find UiGU such that for all 5\J GV 
W (uf , U 1? 5U) + W (uf, U 1; 5U) 
< +W e (f + fi (t) , g + gi (t) , SU) = W a (Ui,5U) (59) 

n 1= n nd (u 1; A M ) + n d (u 1; a m ) 

_ Ui = Uo+Uoi (t) on dV 00 ,u 

where A^ are the Lagrange multipliers computed for the dynamic problem. 
The difference between ( 151?]) and Q5Bj) leads to 

' Find UoiGU such that for all <5U GV 

aw, (n , nf , Uoi, u ,5u) + w (uf, u + u 01 ,8u) 

+ W e (fx (t) , gl (t) , 5U) = W a (Uoi,5U) 

v 7 (60) 

n =n (u , A fcj0 ) , ui = u (u + u i, A M ) 
n 1= n w (u + u 01 ,A M ) + n? 

^ Uoi = Ui (t) on dV 00 ,u 



with 

aw, (no, nf , u i, u ,5u) = w (nf , u + u i,5u) - w (n , u ,5u) (6i) 



knowing that Ui = Uoi- Eq. fl60|) describes the dynamics around a statically 
prestressed configuration for the general case. Observe that all the equations 
are defined in the Lagrangian configuration Vqq = Vo,, free of any external 
prestress effect by definition. The case of moderate rotations is retrieved in- 
troducing in the same equation W- b ^ instead of W, and wf instead of Wf . 
From Eq. (I3T]) and recalling (|36|) . (1521) and (I53I) . Eq. (IBT!) becomes 

-Awf ] (n ,nf,u i,u ,5u) 

= Jv 00 (nf (1) + Uf ( 2p) ) : (V s (5U) + 5 X (U + U i,*U)) dV (62) 

- ; yoo (n« + n^) : (v s (<5U) + 5* (u 0l *u)) ^ 
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8 Moderate rotations and Bernoulli-Navier kinematic assumptions: 
undamped case 



In this section, the stationarity conditions (TTHT) with suited stress constraints 
are used in conjunction with the strain measure (I33p for moderate rotations 
and the so-called Navier kinematic assumptions for beams, in order to obtain 
the corresponding strong form of the dynamic equilibrium equations. This 
analysis will enable a better understanding of the general equations previously 
presented. In particular, a simple way of estimating 77 and p is suggested with 
reference to the example of a clamped-clamped beam. In the non deformed 
configuration Voo = Voi, the beam axis coincides with the cartesian axis X and 
the beam motion is supposed to be limited to the plane X — Y. A quadratic 
Helmholtz energy is adopted, leading to a linear constitutive law depending 
on the tensors D and A: 



II = n 0i + D : E — A (T — T ( 



(63) 



with 



"A 



"A 



E = E< 6 > (U) - £ A k R k =e (U) + X (U) - £ A fc R fc 



(64) 



k=\ 



k=l 



For an isotropic material, one has D 



\v\ ® 1+(1 - 2i/)I]. More- 



over, since A =izf^l, one has D _1 : A (T — T 0i ) = a (T — T 0i ) 1. The Navier 
kinematic assumption reads 



u = [ u -y v',v,o] rj 



(65) 



where u = u (X, t) , v = v (X, t) are the X-and Y-displacement fields; the apex 
' indicates the derivation with respect to X. The strains e (U) and x (U) read: 



s(U) 



u' - Y v" 





X(U) = - 



v /2 
v' 2 




For the prestress n j, we assume: 



n 



0/ 



n 0i ,i 
n 0i , 2 
n 0ii3 



n 



(x) 
Oi 








yriv) I yjjib) 






tt( 2 ) 1 y-rr 
LL 0i t 1 11 ( 



(6) 
0i,z 
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Off-diagonal terms of Iloi may not vanish. Nonetheless, they have no influence 
on the following analysis, since they are associated with zero virtual strain 
components in the virtual work product. The thermal and load fields write 

g x (X,t)-Y g b (X,t) 
9y (X,t) 


T - T 0i = [T x (X) - T 0hX (X)] - Y [7 (X) - 70i (X)] = AT X — Y A7 

The stress constraints usually imposed to retrieve beam equations are Hyy = 
IIzz = 0, formally expressed by the conditions 













Ri : n = R 2 : n = with R x = 


1 


, R2 — 
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Eqs. f j22|) . fl24l) and fl29l written for this case become 

1 = (4 Iv 00 K - Y -"f dV) 1 + Iy 00 3 [a (AT X - Y A 7 )f dv) 1 

= (Wo /o L °° ( A (^') 2 + J K) 2 ) dX) h + (fg /*» (A (AT,.) 2 + J (A7) 2 ) dX) 
V P = (i; fv 00 v' 2 dvf = (^J^A^dxf 

n 0i = o(t 1 ) + o (r] 2 ?) 

< 7] < 1 p> 

where A is the area of the generic beam section; J is the inertia moment and 
Loo is the beam length. Then, the same procedure as in the general case can 
be applied here, in order to determine the strain-rotation domains H a , H& and 
H c . As it is well-known, the Navier-Bernoulli kinematic assumptions entails 
that all shear strains, i. e. the off-diagonal elements of e, are equal to zero. As 
a result, one can easily prove that also the terms of order 3 + p, 2 + p and 
1 + 3p in the strain energy 

riNB 

T = (r, 2 ) + O ( V 1+2p ) + O (r/ 4p ) + E (^) (68) 

become zero (see Figure [8] and compare to Figure [2]). As a result, the strain- 
rotation domains are not the same as in the general illustrated in 



f : 



f x (X,t)-Y f b (X,t) 
fy(X,t) 
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Figure [9j The difference is highlighted by the small region excluded in the 
general case and admitted by the Navier kinematic conditions. Inside this 
region, the pertinence of the Navier assumptions (I65p should be further inves- 
tigated. Other kinematic assumptions, like for instance those of Timoshenko, 
appear to be more sound. 



8. 1 Strong form equations 



Eqs. 



and lead to the Lagrange multipliers 



Ai 



A 2 



n 0i ,2(i-i' 2 )--i/(i+i/)ri w 



+ v (u' - Yv" + ^ 



J2 



\v' 2 - a (AT X - Y A 7 ) (1 + v) 



(69) 



noi, 3 (i-^ 2 )-^(i+^)n i 



+ v(u' -Yv" + W 2 



a (AT X - Y A 7 ) (1 + v) 



from which (see (ED) E XX = E XX = u' - Yv" + ^v' 2 , £ yy = E YY - Ax 



±v /2 - Ai, E zz = E zz - A 2 = -A 2 and 



n X x = noi,, + F n 0i , 6 + £ U - Yv" + V 2 - a (AT X - K A 7 )) 



where n 0i , x = n$ - v (n$ + n$) and n oi , 6 = ng, - z/ (ng, + ng 2 , , 

while the other stress components are zero. Moreover, the virtual works read 



(II, U,<5U) + W e (f , g,5U) = W a (U,5U) 



wi 6) (n, u,5U) = - ; yoo n xx [5u' - ysv + vW] dv 

= wP (n, u,5U) + wf +2p) (n, u,5U) + wf p) (n, u,<su) (70) 

-Wf ) (n, U,5U) = J Lo ° [EA (u' - aATj;)] <fu' c/X + / L °° [£ J (v" - «A 7 )] <5v" 



-W 



(l+2p) 
i 

;(4p) 



n, U,<5U) = f L °° \eA\v' 2 ] Svl dX + J L °° [EA (u' - «ATJ] [v'5v'} dX 



-wr (n, u,5U) = / ( 



Loo 



eaW 2 



[v'SV] dX 



W e (f , g,5U) = / yoo [(/, - F / 6 ) (5u - YSv') + /„ 5v] dV 

+ Uv 00 ,„ [(9* ~ Y 9b) (Su - Y5v') + g y 5v] dA 
+ [RJu + C z 5v' + R y fc}xZo 00 
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W a (U,5U) = J Voo poo [(u - ^ V) (<Ju - W) + v5v] 

= J L °° Pm A n6ndX + J L °° Poo Av5vdX - J L °° Poo Jv"5vdX 



where R x R y and C z are the reaction forces and the reaction moment at 
the boundary <9Voo,«- Observe the second and third term in W a : it can be 
proven that the ratio r\ between the third term (rotational inertia) and the 
second term (translational inertia) reads r\ = O := ^(f?~)> wnere 

h is the beam width. When the squared aspect ratio (h/Loo) 2 is small, r\ 
is small too. Note that this ratio can be easily expressed in terms of r\ and 
7] p when AT X = A 7 = and u' = O (r] 2p ). In this case, Eqs. (1671 -1.2) entail 
r\ = O (r) 2p )+0 (r) 2 ~ 2p ). Hence, under these assumptions r] 2p < ( and rf~ 2p < ( 
suffice to have 7*1 small. The first condition is always fulfilled due to physical 
linearization assumption, while the second one is equivalent to log rf~ 2p < log ( 
and y < x + logC/2 (see Figure [H]) and is satisfied in a large portion of H&. For 
simplicity, the rotational inertia is always omitted hereinafter. The strong form 
equations corresponding to (1761) are derived using the standard procedure: 



AU 0i:X + EA (u' + iv' 2 - aAT x )]' = -p x + p 00 Au 
(EJ (v" - «A 7 ) + m 0i;b y 
_ - (V [AU 0i>x + EA (u ; + \v' 2 - «AT, 



(71) 



-Q 1 +Py~ PooAv 



where p x = Af x and p y = Af y are the horizontal and vertical loading per unit 
beam length, respectively; q = J f b is a couple per unit length. The boundary 
conditions of type <9Voo, CT involving the external forces Px = Ag x , Py = Ag y 
and the external couple M = Jg b at the ends of the beam, are not reported 
for brevity. At the configuration Voo = Vbi, one has u = v = AT X = A7 = 0, 
dVoo = <9Voo,o- with zero loads, entailing (AHoi tX ) = and (JTIoi^) = 0. This 
implies that AIlo^x is constant along the length of the beam. Since zero force 
is applied on dVoo, i.e. Px = Ag x (0) = = ATIoj^, this constant is equal to 
zero and the same holds for Jlloi^. This means that the initial self- equilibrated 
stress is zero and the equilibrium equations become 



¥ 2 



- (EA (u' - aAT x 
(EJv" - aEJA-f)" 
- (y [EA (u' + |v' 2 - a AT 



p x - poo Ail 

= -q' +p y - pooAv 



(72) 
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This is the general expression of the beam equation with temperature field. 
The corresponding expression of the strain energy (see (131]) and (1361) ) becomes 



^(2) = jLoo EA fl u , _ aAT ^ u , rfx + jLoo £ j Q y // _ a A ^ v „ dx 
jr(i+2 P ) = jioo EA (l u > _ aAT ^ 1 y /2 rfx + 1 jioo EA Q v ' 2 ) u' dX 



S.S T/ie geometric interpretation of rj andrf 



In this Section, an interpretation of r\ and ?7 P in terms of suitable deflection 
and shape ratios, and as functions of the temperature field is provided for 
the case of a homogeneous beam. A first example concerns a beam with very 
small bending stiffness, i.e. J/ A ~ 0. A static vertical load F is applied at 
the midspan, where it induces a transversal displacement v max . Moreover, ii = 
v = p x =Py = q = and an axial temperature field is introduced. Then, Eq. 
(1721) becomes 

EA (u' - aAT x + ±v' 2 ) = R x 
R x v" = 



where R x is the constant horizontal reaction at X = Lqq. Boundary conditions 
write u (0) = 0, u (Lqq) = u > 0, v (0) = v (L 00 ) = and v" (0) = v" (L 00 ) = 0. 
Since v' is piecewise constant, with a discontinuity at the midspan, integration 
of the first equation in ( 1731) yields 

Kx - U + -v' 2 - — t aAT x dX =— + -v' 2 - aAT x 



EA Lqq 2 Lqo J -^oo 



where AT X is the averaged temperature variation. It follows u' = aAT x — 
aAT x + i.e. u' = when the axial temperature field is constant, even if 
non-zero. In this case AT X = AT X = const, and by using (1671) . one obtains 



V = Ve + VAT = ^- + V3a\AT x \ and rf = |v'| = 2^^- (75) 

-^00 -^00 



which provide a simple interpretation of rj and t] p in terms of temperature 
difference and of ratios between the maximum displacements and the beam 
length. Hence, the strain-rotation domains of Figure which depend on rj and 
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r] p , can also be interpreted using these ratios. For instance, consider the case of 
zero temperature field and given u value, such that r\ = r] £ = u/L 00 = 10~ 8 : the 
strain-rotation points corresponding to this situation and for different values 
of |v max | are depicted in Figure [U]-b: they have the same y- value (constant 77) 
and different x- values. The larger |v max |, the larger rf: then, according to the 
value of |v max | , the point representing the structural state may belong to any 
of the sets H a , Elf, or H c and the relevant equilibrium equation is different in 
each case. Since the strain-rotation domains depend on (, Figure [9] refers to 
the case C = 0.01. Observe in addition that A a : = -k=" !2AT - = Q ^ T ^ gives an 

estimate of the relative importance of the thermal and mechanical strains. 

Let us now consider a homogeneous beam with a distributed vertical static 
load p y , a generic temperature field and with u = v = p x = q = 0. The same 
structure has been studied in the numerical examples of Section 5.3. Eq. (1721) 
becomes 

I EA (u' - aAT x + Iv' 2 ) = R x 
[EJ(V" -oA 7 ")-V =Py 



For the X-direction, the boundary conditions (b.c.) are u (0) = u (Loo) = 
and in the vertical direction one has v (0) = v (L Q0 ) = v' (0) = v' (Loo) = 0. 
Integrating the first equation with the b.c. at X = and X = Lqo leads to 



1 ^oo 1 ^oo 1 

- / aAT x dX + — ( - v' 2 dX (77) 
'oo J -t^oo J <Z 



EA Lqo j -^00 



It follows, according to (1751-1) 



1 L T\ 1 



u' - aAT x = — / -V 2 dX - -V 2 - aAT x (7S 

iy J * ^ 







where AT X has the same definition as in the previous example. Assume that 
AT X and A7 are constant and substitute ( 1781) into the definition (1671) of 77. 
Hence 

V = Ve+ VAT, Ve = J? 00 (a - + J (v*) 2 ) ctt) " 

VAT =(3a 2 ((AT x ) 2 + i(A 1 ) 2 )Y 



In order to have a better understanding of the geometrical meaning of 77, the 
solution v (X), depending onp^, AT X and A7 should be analytically expressed. 
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However, this is not a simple task in general. Hence, accounting for the b.c, 
we assume here that the deformed shape is approximately co-sinusoidal: 



y(X) 



cos 



/2ttX 
[~L~ 



00 



(79) 



Hence, by using the definitions ([57]) . one obtains 



p TT |vmax| 

' V2 L 00 ' 



V = Ve + Vat 



t 2 [Vmaxj / vg 



Ve=t^ 1 T T\lTt + Q4C L^ 



h 2 



Vat = VSaJAT^ + c (A 7 /i) z 



where c = -ttj, with ft the beam width. When the transversal displacement is 
different from zero, the ratio between the thermal and mechanical contribu- 
tions in rj is equal to 



A a :- 



Vat 



a 



c(A 7 /i) i 



The interest of Eq. ( IHUl) is that it gives a geometrical interpretation for the 
case of beams of the quantities rj and rf defined in ( 122]) for the general case and 
in ( 1671) when the Navier-Bernoulli kinematic assumptions are adopted. They 
can be easily related to geometrical ratios involving the maximum deflection, 
the width and the length of the beam. These geometrical ratios are known 
to be important for beam analysis, but they are clearly related here to the 
tensorial quantities of a full 3D formulation of the structural problem. When 
both maximum deflection and geometry of the beam are known (or estimated), 
it is easy to find the corresponding point in strain-rotation domains and to 
establish how many terms need to be taken into account for the computation 
of the solution. 



In order to find the explicit expression of the strain energy, note that AT X 
constant entails from Eq. 



u 



1 

-^00 



^oo . 



r* d X - -v' 2 



J2 



(82) 
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and Su' = f] p 5r] p — v'Sv'. Hence, using (1791) and the expression of r] p given in 
( 1B01 . the strain energy contributions (173"]) read 



h 

. Loo 



2 /.. \ 2 

max 



JW=EAL M (fife) +vr 4 c( 
-77(1+2*) = eALqq (-fj (^) 4 - ^aAT,. (W 2 



The ratio between the energy terms of order 4p and 2 reads 



•^ (4P) _ o W = 3 



:F(2) 64c +(^) 2 8 

and it only depends on the ratio v max //i. The virtual works of internal and 
external forces (see (1701) ) read 

-W t (II, U,5U) = EALqq U (^) 3 + 2vr 4 c (4) 2 fe ) - 2f aAT. ( 



-^00 



W e (f,g,5U)=^°> 2/ 



l-cosf^l dX = ^Sv- 



\ Loo 



2 



1 max 



with p y supposed constant, and by the virtual work principle, one obtains 

-aAT x ) + (Y^) 3 = hbl (83) 
2 I L 00 8 VLqo/ 2EA V ; 




which provides a simple nonlinear relationship between v max and p y . Eq 
does not depend on A7, since this quantity is assumed constant on the beam 
length. Using Eq. ([83]) with AT X = and supposing h/L 00 = 0.01/0.5 = 0.02, 
c = 1/12, A = lcm 2 and p y = Afyo daN/cm like in Section 5.3, one obtains 
the v max values collected in the first column of Table [7J Moreover, by means of 
( 1501 . it is easy to compute r\ = r/ e and rf . These values are good estimations 
of the corresponding quantities computed from the numerical tests (Table [3]). 
The difference between 77 estimated as in (IHUl) and in the numerical analysis 
is due to the shear stress as well as the Hyy component. On the other hand, 
r] xx issued by the numerical analysis and the analytical estimation of 77 = rj e 
according to (ISO"]) are very close. The corresponding co-ordinates in the strain- 
rotation domain are reported in the fifth and the sixth column of Table [7J 
(see also Figure [9]). The influence of a thermal field with a = 10~ 5 °C _1 , 
AT X = 20° C and hA^ = 10° C on the strain can be estimated as follows: Eq. 
( 1S31 with the true temperature field is used to compute v max , then r\ e and 77 at 
are separately evaluated, as well as their ratio A a (see the seventh column of 
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Table [8]). A more useful comparison may be done by the ratio between v max 
of this last case and v max for the case AT X = at the same level of p y (see the 
last column of Table [SJ). This ratio correctly neglects the influence of A7. 



8.3 Equilibrium around a prestressed configuration: strong form equations 



By subtracting Eq. (JTzj) written at the state V\ from the same equation at the 
state Vq, one obtains 



- [EA (u' 01 + + v^vyj = p%, x - Poo^hoi 
(EJV \)" - (Voi^no,)' 

'oi + v 'o) [ EA ( u 'oi + ¥'m + v i v o)] ) = ~Qi + Pi,v ~ Poo^v 01 



where uo, v are the axial and transversal displacements characterizing the 
statically prestressed configuration, measured between Voo and V ; T 0>x — T 0itX 
is the axial temperature field; H 0iX = EA (u' — a (T ^ x — T 0itX ) + |vq 2 ^) is the 
static axial prestress; the unknowns u i = Ui — u , v 01 = v 1 — v are the 
displacements between the dynamic configuration V\ and the static one. The 
first and the second term of the bending equation (|841 -2) are linear. The second 
one is related to the static configuration Vq (uo, v and To }X )' it is the effect of 
the static prestress due to external static loads, i.e. it is not self-equilibrated 
and is associated with a body configuration Vq different from Vqi = Vqq. The 
bending equation in (1511) contains a term coupling u' 01 with the static and the 
dynamic rotations v' and v' Q1 . This term seems to be of paramount importance 
since experimental investigations flTreyssedd . 120061 ) prove that pre-bending 



effect v' play a crucial role, and not only prestress. When u' 01 is different from 
zero, this term cannot be neglected in front of the static prestress contribution 
associated with U 0x . This situation occurs when a pulsating axial loading 
Px (t) is applied at one end of the beam, inducing parametric resonance. In 
this case, the first equation in becomes 



EA 



u 'oi + ^ v 0l + v r 



01 v 



-Px(t)+ J p 00 AuoidX 



where p\ x = by ass umption . The i nfluence of the axial inertia on the solution 
is discussed, e.g., by iRibeirol (]200ll ). 
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9 Moderate rotations and Bernoulli-Navier kinematic assumptions: 
viscous damping case 

The stress when linear viscous damping occurs is given by (see (j47p . (j49]) and 

my- 

U=U nd +U d U nd = n 0i + D : (E- E^i AfeRfe) -A (T - T oi ) 

where the Lagrange multipliers can be computed imposing the constraints 
R fc : U nd = and R fc : n d = k = 1, 2 , where Ri and R 2 are given in (1551) . 
Since the two parts of the stress must be separately equal to zero, the Lagrange 
multipliers A 1; A 2 do not change with respect to the non-dissipative case, and 
they are given in Eq. (|69|) . Hence, their time derivatives can be computed 
under the assumption, previously discussed, that the time variations of the 
temperature field are very small compared to those of displacements. Hence, 
Ax = us + v'v' and A 2 = us, with s = u' — Yv" + vV. As a result, one has 

s 
-us - Ai 
-us - A 2 

Let us assume that damping is proportional to the stiffness and to the mass, 
viz. F=a 5 p 00 I+/%D =a ( pool+fc (Al ® 1+2/iI) = A d l ® l+2// d I, with A d = 
/^A and 2/i d = «g p o + 2/i/?^. Then, the constraints on the dissipative stress 
read 

( F ■■§) YY = Poo (-m - Ai) + P Hl+u) l_ 2v) [(1 - i/) (-Ax) + u (-A 2 )] = 
( F --%zz = a Z (~ PS - X ^ + [" ("AO + (1 " ") (-Aa)] = 

It entails A, = A 2 = _^^ a , U d YY = U d zz = and U* xx = 
(F :f ) xx = c ?S with 

= (&£) 2 + q € pop^g (2 - z/) + (1 + i/) (1 - 2z/) (q g p 00 ) 2 
° 5 + (1 + i/) (1 - 2u) ^ poo 

This expression relates the "structural" viscous damping coefficient eg and 
the "material" parameters a^,^, poo, E and u. The virtual work of dissipative 



dE 

~dt 



dE 

~dt 



fe=l 
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forces becomes (see Eq. (J56 



W? = Wi (ll d , U,5U) = - J U d xx (5u' - Y5v" + v'5V) dV 

Voo 



and the strong form equations read 

- EA (u' — a (T x — T oitX ) + |v' 2 + c 5 (u' + vV 
(EJ ( v "-a( 7 - 7oi ) + Ce v"))" 



Pa - Poo-4u 



55) 



£A u' - a (T x - T 0hX ) + ±v /2 + c c (u' + vV) = -q> + Py - Poo Av 



When F=/3^D, i.e. damping proportional to the stiffness, one has Ai = A2 = 
and = (3^E. Moreover, around a static configuration with prestress Hq x = 
E (u' — a (T 0x - T QitX ) + |v' 2 ) , one has 



- EA (u' 01 + I v' 2 + v(,v' 01 + c c (u' 01 + (v' + v' c 
(EJ (v^ + c^))"-^^)' 

A 01 + V ) [jM (<! + iv i + V' V' 01 + C ? (i^ + K + V^ 

= ~q[ + Piy - Poo-4voi 



Pix - Poo^hoi 



01 J v 01 



(86) 



10 Conclusions 



Based on a Hu-Washizu functional accounting for stress constraints, an orig- 
inal derivation of the dynamic equations for thin prestressed and prestrained 
structures from continuum mechanics has been presented. The assumption of 
small strains and moderate rotations, as well as the physical linearization, have 
been formalized by the notion of strain-rotation domains associated with new 
strain and rotation global measures. The corresponding approximated expres- 
sions for the strain energy have been justified by 2D non-linear finite element 
investigations. The presence of a generic mechanical or thermal prestress and 
physically linear viscous damping has been discussed in detail. In particular, 
damping is introduced by the original definition of a pseudo-potential with 
stress constraints. Coupled beam equations governing traction and bending 
for small strains and moderate rotations have been derived by using the pro- 
posed general procedure. 
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11 Appendix: Lagrange multipliers 

In this Section, we show how to compute the Lagrange multipliers. The stress 
constraints write 



or K A = B, where A= (A;) i=1>nA , K = (%) ii=lj7lA , B = (Bi) i=1>nA and 
k i:j = R 4 : D : R, , B { = R: (U 0l + D : E (U) - AAT) 




(87) 



with k 



1,71a- Eq. flHTl) is equivalent to 



£ R fe : D : A; = R k : (II 0l + D : E (U) - AAT) 



k=i 
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By virtue of the symmetry and positivity of D it is easy to prove that K is 
symmetric positive definite and therefore invertible, provided that tensors 
are linearly independent. Hence, the Lagrange multipliers read A = K -1 • B. 
Then, the assumptions §M) and flS} prove that A = A (1) + + 

Ay-+p) 4- A*- 2 -*. In particular, the constraints fl66|) leads to 



K 



A + 2/x A 
A A + 2/i 



and 



Ai 
A 2 



(A + 2nY - A 2 




{U i + D : E (U) - AAT) yy 
(n , + D : E (U) - AAT) ZZ 



For the dissipative stress, one has analogous constraints 



'dE(U) ^dAi 



II A 



Then A = K 1 B = A (1) + A (2p) + A (1+p) + A (2) , where A = (A;) J=1 nA and 



K = 

B = , 5, = R: F : (™ - ^R) 



K is symmetric positive definite. For a beam 
/ 



K 



Ad + 2/id Xd 
\ A d A d + 2// d 



and 



Ai 
A 2 



(A d + 2fji a 



\ 2 



Ad + 2/id —Ad 
~Ad Ad + 2/i d 



rp . / rfE(U) Y^ n A rfA; I > 
1 ' 1 _ It ^1=1 ~dT t ^ 1 

* • 1 — It L>i=i ~dt^j 
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Fig. 2. Orders of the strain energy terms. The dominant terms (2, l + 2p and 4p) are 
highlighted. The line corresponding to the order 6p represents the most important 
energy term associated with the non-linear constitutive law ([8]). 
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Fig. 3. (a) Strain-rotation domains for a physically linear material, (b) Strain-rota- 
tion domains under the assumption of physical linearization. The sets W b and 
are merged into the set The "strain-rotation points" represent the deformed 
configurations of the structure of the numerical example discussed in Section 5.3 . 
In both Figures (a) and (b), the dominant terms of the strain energy are indicated. 
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10.01 
10.05 
1 0.5 



Fig. 4. Volume vertical force fyo = — 1 daN/cm 3 . 
* (E) [daNcm cm,- 3 }, with E = E (U). 



Map of the strain energy density 
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IsoValue 

■ o 

■ 0.01 

■ 0.05 

■ 0.5 



Fig. 5. Volume vertical force fyo '- 
energy density (E) [daNcm cm~ 



-1 daN/cm 3 . Map of the approximated strain 
with E = E^ (U)). 
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IsoValue 

■ -0.00065 

■ -0.0003 

■ o 

■ 0.0003 

■ 0. 00065 



Fig. 6. Volume vertical force fyo = — 1 daN/cm?. Map of the relative difference be- 

, , , ■ , <I'(E(U))-*(E( i, )(U)) 

tween the exact and the approximated strain energy density, viz. #(e(u)) ■ 
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IsoValue 

■ -0.00065 

■ -0.0003 

■ o 

■ 0.0003 

■ 0. 00065 



Fig. 7. Surface vertical force gyo = — 1 daN/cm 2 . Map of the relative difference be- 

, , , ■ , <I'(E(U))-*(E( i, )(U)) 

tween the exact and the approximated strain energy density, viz. WTeTuY) ' 
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0.2 0.4 0.6 0.8 1 

P 

Fig. 8. Navier kinematic assumptions and orders of the strain energy terms. The 
dominant terms (2, 1 + 2p and 4p) are highlighted. The line corresponding to the 
order 6p represents the most important energy term associated with the non-linear 
constitutive law (j£J) . Note that the terms of order 3 + p, 2 + p and 1 + 3p no longer 
appear (compare with Figure [2]). 



43 




Fig. 9. Physical linearization and Navier kinematic assumptions, (a) Strain rota- 
tion domains, (b) The strain-rotation domains, where M' b and W£ are merged. The 
"strain-rotation points" represent the deformed configurations of the structures of 
the examples discussed in Section 8.2 . 
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Ui = U + Uoi 
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E(U ) 


E(Ui) 
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<5E (0,<JU) = 5e 


5E(XJ ,5XJ) 


<5E(Ui,<5U) 
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n () = n 0i + D : E(U ) 


n, H +D :E(Ui) 




-A(To-Toi) 


-A(To-Toi) 


f 





fo 


fo + fiW 


g 





go 


go + gi(*) 


u 





Uo 


U + Ui(t) 


u - 








tii = Uoi 


Table 1 









Displacements, strains, temperatures and stresses associated to each structural con- 
figuration. Uoi = Uoi = in the analyses of this paper, since no geometric imper- 
fection is considered. 
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(n) \ (E) 
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T 

r • r 


(e -r) 




e, aAT 


2 


1 + 2p 


2+p 


3 


T 

r • r 


1 + 2p 


4p 


1 + 3p 


2 + 2p 


T 

e 1 ■ r 

T 


2 + p 


l + 3p 


2 + 2p 


3+p 


T 

£ ■ £ 


3 


2 + 2p 


3 + p 


4 


£ ■ £, (aAT) 2 


3 


2 + 2p 


3 + p 


4 


(r ■ rj • (r • rj 


1 + 4p 


6p 


1 + 5p 


2 + 4p 


(f . r ) . ( e T . r ) 
(r T • e) • (r T • e) 


3 + 2p 


2 + 4p 


3 + 3p 


4 + 2p 


i. ■ e) ■ ( • e) 


5 


4 + 2p 


5 + p 


6 


(e) • (r T • r) , aAT (r T • r) 


2 + 2p 


1 + 4p 


2 + 3p 


3 + 2p 


e • (e T • r) , aAT (e T • r) 
£ ■ (r T ■ e) , aAT (r T ■ e) 


3 + p 


2 + 3p 


3 + 2p 


4 + p 


£ ■ £ T ■ £, aAT £ T ■ £ 


4 


3 + 2p 


4 + p 


5 


(r T ■ r) • (£ T ■ r) 
(r T • r) • (r T • e) 


2 + 3p 


1 + bp 


2 + 4p 


3 + 3p 


(r T • r) • ( £ T ■ e) 


3 + 2p 


2 + 4p 


3 + 3p 


4 + 2p 


(£ T • r) • (£ T ■ £) 

{r T • e) • (e T • e) 


4 + p 


3 + 3p 


4 + 2p 


5 +p 



Table 2 



Orders of the terms appearing in the strain energy expressions (|30j) -()32|) and (|421 2). 
representing the physically linear and non-linear case, respectively. 
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\fv,o\ 
[daN/cm 3 ] 


I Vmax 1 

[cm] 


max | £ xx | 


Vxx 




V 




if 




X 


y 


P 


0.1 


0.0092 


6.46e - 5 


1.52e - 


5 


1.59e 


- 5 


4.05e - 


4 


-3.4 


-4.8 


0.71 


0.5 


0.0459 


3.20e - 4 


7.58e - 


5 


7.93e 


- 5 


2.02e - 


3 


-2.7 


-4.1 


0.66 


1 


0.0914 


6.34e - 4 


1.51e - 


4 


1.58e 


-4 


4.03e - 


3 


-2.39 


-3.8 


0.63 


3 


0.263 


1.85e - 3 


4.37e - 


4 


4.66e 


-4 


1.157e - 


- 2 


-1.94 


-3.33 


0.58 


6 


0.4745 


3.58e - 3 


8.03e - 


4 


8.96e 


-4 


2.09e - 


2 


-1.68 


-3.05 


0.55 



Table 3 



Clamped-clamped beam (2D plain stress analysis). Numerical results for different 
values of the vertical volume load fyo- 
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\fv,o\ 




P2 = V Ap - 2 


T 




T 


[daN/cm 3 ] 






[daN cm] 


[daN cm] 




0.1 


0.0103 


0.00011 


0.00613653 


0.00613653 


~ 


0.5 


0.0515 


0.00265 


0.153092 


0.153091 


6.5e — 6 


1 


0.1025 


0.01052 


0.608433 


0.60842 


2.1e - 5 


3 


0.2873 


0.0825 


5.15335 


5.15244 


1.8e-4 


6 


0.4875 


0.238 


17.8727 


17.8626 


5.7e-4 



Table 4 



Clamped-clamped beam (2D plain stress analysis). Energies for different values of 
the vertical volume load fyo- 
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\9Y,o\ 
[daN/cm 2 ] 


I Vmax | 

[cm] 


max \ e xx 1 


Vxx 




V 




rf 




X 


y 


P 


0.1 


0.0092 


6.46e - 5 


1.52e - 


5 


1.59e - 


5 


4.052e 


-4 


-3.39 


-4.80 


0.71 


0.5 


0.0459 


3.20e - 4 


7.58e - 


5 


7.93e - 


5 


2.02e - 


- 3 


-2.69 


-4.10 


0.66 


1 


0.09142 


6.34e - 4 


1.51e - 


4 


1.58e - 


4 


4.029e 


- 3 


-2.39 


-3.80 


0.63 


3 


0.263 


1.85e - 3 


4.37e - 


4 


4.67e - 


4 


1.158e 


- 2 


-1.94 


-3.33 


0.58 


6 


0.475 


3.58e - 3 


8.04e - 


4 


8.98e - 


4 


2.092e 


- 2 


-1.68 


-3.05 


0.55 



Table 5 



Clamped-clamped beam (2D plain stress analysis). Numerical results for different 
values of the vertical surface load gyo- 
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\9Yfl\ 


Pl = if?- 1 


P2 = V Ap - 2 


T 




T 




[daN/cm 2 ] 






[daN cm] 


[daN cm] 






0.1 


0.01035 


0.00011 


0.0061368 


0.0061368 


~ 




0.5 


0.0515 


0.00265 


0.153124 


0.153123 


6. 53e - 


6 


1 


0.1025 


0.0105 


0.608689 


0.608676 


2. 14e - 


5 


3 


0.2870 


0.0824 


5.15922 


5.15831 


1.76e - 


4 


6 


0.4876 


0.2377 


17.9055 


17.8953 


5. 70e - 


4 


Table 6 



Clamped-clamped beam (2D plain stress analysis). Energies for different values of 
the vertical surface load gyo- 
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\Py\ 
[daN/cm] 


V max ,AT=0 

[cm] 


V 








a; 


y 


P 


0.1 


0.009165 


1.477 e 


-5 


4. 072e - 


4 


-3.39 


-4.83 


0.702 


0.5 


0.04576 


7.377 e 


- 5 


2.033e - 


3 


-2.69 


-4. 13 


0.651 


1 


0.0911 


1.469e 


-4 


4.047e - 


3 


-2.39 


-3.83 


0.624 


3 


0.2616 


4. 243 e 


-4 


1.162e- 


2 


-1.935 


-3.37 


0.574 


6 


0.4714 


7. 754 e 


-4 


2.094e - 


2 


-1.68 


-3.11 


0.54 


Table 7 



Clamped-clamped beam with co-sinusoidal deformed shape. Numerical results for 
different values of the vertical load p y . 
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\Py\ 
[daN/ cm] 


I Vmax | 

[cm] 








VAT 




'/ 




rf 






| v max | 
| v max,AT=o| 


0.1 


0.010808 


1. 


. 742e 


- 5 


6. 083e - 


4 


6.26 0e 


-4 


4. 802e 


-4 


34.9 


1. 180 


0.5 


0.053906 


8. 


690 e 


- 5 


6. 083e - 


4 


6.952e 


-4 


2.395e 


- 3 


7.0 


1.178 


1 


0.1070 


1. 


726 e 


-4 


6. 083e - 


4 


7. 809 e 


-4 


4. 754e 


-3 


3.52 


1.174 


3 


0.3003 


4. 


8801e 


-4 


6. 083e - 


4 


1.096e 


- 3 


1.334e 


- 2 


1.25 


1.148 


6 


0.52242 


8. 


632 e 


-4 


6. 083e - 


4 


1.472e 


-3 


2.321e 


- 2 


0.705 


1.108 



Table 8 

Clamped-clamped beam with co-sinusoidal deformed shape. The vertical load is the 
same as in Table [3l Moreover, a thermal field with AT X = 20° C and /1A7 = 10°C 



has been added. 
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